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TECHNICAL NOTE 3331 


ANALYSIS OF LAMINAR FORCED-CONVECTION HEAT TRANSFER IN ENTRANCE 
REGION OF FLAT HECTANGDLAR DUCTS 
By E. M. Sparrow 


SUMMARY 


The simultaneous development of temperature and velocity profiles 
in the entrance region of a flat rectangular duct is studied. The flow 
is assumed to he laminar with negligible dissipation. All fluid proper- 
ties are taken to he constant , and the wall temperatures are uniform. 
Thermal and velocity boundary layers are calculated using the Karm^n- 
Pohlhausen method. Nusselt numbers are reported for Prandtl numbers in 
the range 0.01 to 50. 


When the Nusselt numbers are plotted against Graetz number 
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there is a separate curve for each Prandtl number. In contrast, the 
Nusselt numbers for the case of a parabolic velocity profile throughout 
can be plotted against Graetz number as a single curve which is appli- 
cable for all Prandtl numbers . Beyond the position where the boundary 
layer analysis can no longer be used, curves have been faired to connect 
the Nusselt number results of the present analysis -with those of Norris 
and Streid, who assumed a parabolic velocity profile throughout the 
length of the duct. 


Calculations were also made for a duct having one wall at uniform 
temperature and the other wall insulated. Nusselt numbers are reported 
for the range of Prandtl numbers from 0.01 to 50. 

The results of the present analysis are compared with those for a 
single flat plate. For purposes of comparison with the results of Norris 
and Streid, the boundary layer analysis was also used to calculate Nusselt 
numbers for the case of an unchanging parabolic velocity profile. 


INTRODUCTION 

Until recently, analytical studies of laminar forced-convection 
flow in tubes and ducts have been concerned primarily with three groups 
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of problems. In. the first group, a fluid with a velocity profile that 
is already fully developed enters a section of pipe having a wall tempera- 
ture different from the temperature of the entering fluid. The initial 
analysis for the development of the temperature profile was made by Graetz 
In 1883 for a round tube. In 1910, Graetz' s solution was obtained inde- 
pendently by Eusselt . Since then, a number of investigators have recal- 
culated and extended Graetz's results; for example, Etrew, Jakob, Lee, and 
Groeber . A different approach to Graetz ' b problem was made by Ldveque 
In an effort to obtain better numerical values in the entrance region of 
the tube section, where Graetz's results are inaccurate. A presentation ^ 

of the aforementioned investigations is given in reference 1 and refer- to 

ence 2. The problem analogous to that of Graetz for flow between two 
parallel planes (flat rectangular duct) has been treated by Korris and 
Streid (ref. 3) and by Prins, Mulden, and Schenk (ref. 4). 

The second class of problems is concerned with the development of 
the velocity profile in the entrance region of a pipe with no heat trans- 
fer. The fluid is assumed to enter the pipe with a uniform velocity. 

In the course of its flow through the pipe, the fluid is retarded by 
friction at the pipe wall, so that the velocity profile is distorted and 
finally becomes parabolic far from the entrance. For the round pipe, 
the problem has been studied by Schiller, Atkinson and Goldstein, 

Boussinesq (all reported in ref. 5, pp. 301-308), and Lahgharr (ref. 6). 

The first three used the boundary layer concept in at least part of the 
analysis. The analogous problem for the flat duct was treated by 
Schlichting (ref. 7) and by Schiller (ref. 5, p. 309), both of whom 
assumed the existence of boundary layers. 

In the third group are so-called fully developed problems in which 
the heat-transfer and friction parameters do not change along the length 
of the duct. The results are applied to long ducts. Examples of fully 
developed solutions may be found in the work of Clark and Kays (ref. 8), 
where the flow and heat transfer in noncircular ducts is studied for the 
cases of uniform wall temperature and uniform heat flux. 

Only recently has some attention been given to the case of simul- 
taneous development of velocity and temperature profiles in the entrance 
region of a pipe. Solutions for air flowing in round tubes have been 
obtained by Kays (ref. 9) and by Toong (ref. 10). With the assumption 
that air is an incompressible, constant-property fluid, a direct numeri- 
cal integration of a simplified energy equation by the use of finite 
differences is obtained in reference 9. The velocity values required 
in this integration were taken from the results of reference 6. The 
analysis of reference 10, which is for high-speed flow, assumes that the r 

air is compressible while the other properties are constant. The boun- 
dary layer differential equations are transformed and then solved on a 
differential analyzer for specific values of inlet Mach number and wall- 
to-inlet temperature ratio. 
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In the present analysis, simultaneous development of the velocity 
and temperature profiles in the entrance region of a flat rectangular 
duct is studied. Incompressible, constant -property fluids having Prandtl 
numbers in the range 0.01 to 50 are considered. The presence of veloc- 
ity and thermal boundary layers is assumed, and these boundary layers 
are taken to have definite thicknesses, which are calculated by the 
Earman-Pohlhausen method,. Nusselt numbers are computed for the length 
of duct over which the boundary layer analysis is made. Beyond the po- 
sition where the boundary layer treatment can no longer be used, curves 
have been faired to connect the Nusselt number results of the present 
analysis with those of reference 3, in which a parabolic velocity pro- 
file is assumed throughout the entire duct length. 

In addition to the symmetrical situation in which both duct walls 
have the same uniform temperature, a second case is treated, in which one 
wall is Insulated and the opposite wall is held at a uniform temperature. 

v 

The Nusselt numbers obtained in the present analysis should also be 
applicable to annuli in which the curvature effects are small (i.e., where 
the ratio of the diameters of the concentric cylinders forming the walls 
of the annulus Is close to 1), 

So that a comparison, could be made with the results of reference 3, 
the boundary layer analysis was used to study the development of the 
temperature profile associated with a parabolic velocity profile through- 
out. The development of the temperature profile was also calculated for 
the case of a uniform velocity profile throughout ( slug flow) .' 

ASSUMPTIONS , 

The physical model selected for a study of the development of the 
velocity and temperature profiles in the entrance region of a flat rec- 
tangular duct is the following: 

(1) The flow is laminar. 

(2) All fluid properties are constant. 

(3) The velocity and temperature profiles are uniform across the 
entrance section. 

(4) There exists a velocity boundary layer of definite thickness 5. 
As on a flat plate, the action of. viscocity in the entrance region 

of the duct is confined primarily to che fluid layers near the wall, 
although in actuality small effects exist everywhere in the flow. It 
will be assumed in this analysis that viscosity plays a role only in a 
definite region adjacent to the wall called the velocity boundary 
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layer S. The fluid outside the velocity houndary layer is called the 
core. Inasmuch as the fluid adjacent to the -wall is retarded by fric- 
tion, it is necessary that the flow in the core he accelerated in order 
that the same mass pass through every cross section of the duct. The 
velocity profile is taken to he flat throughout the core. The velocity 
boundary layer thickness is assumed zero at the entrance section. The 
boundary layer grows in thickness along the length of the duct until it- 
reaches the center line, where it meets the boundary layer from the 
other •(■rail of the duct. Figure 1 shows diagrammatically the development 
of the velocity boundary layer. ^ 

to 

(5) Two temperature conditions at the duct walls are to be 
considered : 

(a) Both duct .walls have the same uniform temperature, which 
differs from the temperature of the fluid entering the duct; this 
will be called the case of heat transfer at both walls. 

(b) One of the duct walls is maintained at a uniform tempera- 
ture which differs from the entering fluid temperature while the 
other wall is insulated; this will be called the case of heat trans- 
fer at one wall. 

(6) The viscous dissipation and work of compression are negligible 
compared with the heat conduction. 

(7) There exists a thermal boundary layer of definite thickness A. 

Since the temperature of the fluid entering the duct differs from the 
temperature of the duct walls, there will be a heat transfer to (or 
from) the fluid. The principal effects of this heat transfer will be 
felt by the fluid layers close to the wall, although small effects will 
exist everywhere. The assumption will be made that the effects of the 
heat transfer play a role only in a definite region adjacent 'to the wall 
called the thermal boundary layer A. The fluid outside the thermal 
boundary layer will be uninfluenced by the heat transfer and will there- 
fore have a uniform temperature identical to the value at the entrance 
of the duct. The thermal boundary layer thickness will be assumed zero 
at the entrance section. For the symmetrical case where both duct walls 
have the same uniform temperature, thermal boundary layers will grow in 
thickness along the length of the duct in an identical manner for both 
■trails. These boundary layers will meet at the center line, where A = a. 

For the case where one duct wall is maintained at a uniform temperature 
and the opposite wall is insulated, a thermal boundary layer will develop 
only at the wall at which the heat transfer occurs. The boundary layer 
will grow until it reaches the opposite wall, where A = 2a. In general, 
the velocity and thermal boundary layers may be of different thicknesses. 
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(8) The flow is two-dimensional. 

(9) The flow is steady. 


BASIC EQUATIONS 


The laws of conservation of momentum, energy, and mass for steady 
laminar boundary layer flow with constant properties and negligible dis- 
sipation and work of compression are as follows: 


u 






afu 

Sy 2 



(la) 

(lb) 


u 





du dv 

s + ^ = 0 


( 2 ) 

( 3 ) 


(The symbols used herein are defined in appendix A.) Eor the core flow, 
which has been assumed nonviscous and uniform across the section, the 
momentum equation reduces to 


U 


1 



( 4 ) 


The pressure term is eliminated from equation (la) by means of equations 
(4) and (lb). Then, using the assumption that definite values A and 
5 are associated with the thermal and velocity boundary layers, respec- 
tively, equations (la), (2), and (3) are integrated across the section to 
yield 



This integration may be found in reference 5, page 136 and pages 614-615. 



6 


MCA TH 3331 


T (u - U )dy + a(U, - U) = 0 (7) 

u 0 

These equations require only that the conservation laws he satisfied for 
the cross section as a whole . This is a much less stringent requirement 
tjian that imposed by equations (l), (2), and (3), which ask that the con- 
servation laws he satisfied at every point in the cross section. 

According to the Karmen -Pohlhausen method, the velocity and tempera- 
ture are approximated hy polynomials in y having coefficients that are 
functions of x. The coefficients are determined hy satisfying boundary 
conditions at the duct wall and at the edge of the boundary layer and 
hy using equations (5), (6), and (7). 

A simple and convenient expression for the velocity profile in the 
boundary layer which gives reasonably good agreement with Schlichting's 
more exact results is the following one first used by Schiller: 


0< y< 6 (8) 


- 4 ©-® 2 ] 


This equation satisfies the conditions 

u = 0 at y = 0 


u 


U-, 



at y = 5 


O) 


but does not satisfy the condition that at the wall (y — 0), 


d 2 u ^ 1 dp 
dy 2 H 


This requirement arises from the evaluation of equation (la) at the wall, 
where u = v = 0. 


For the temperature profile in the boundary layer, the following 
relation is chosen: 


t ~ S = 3(Z\ 

H ~ \ 2 W 



0 < y < A 


( 10 ) 


The boundary conditions satisfied by this expression sire 


t 



dy 2 



it y = 0 


( 11 ) 
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The condition, 

3y 2 

(2) at the wall. 


0 



at y = 0 arises from the evaluation of 


(ID 

Cont . 
equation 


It is seen from equations (8) and (10) that 8, Hi, and A remain 
to he determined as functions of x. .The assumption of constant fluid 
properties means that the solution for the velocity (i.e., 8 and U^) is 
independent of the heat transfer. (The velocity solution therefore will 
apply for "both the case of heat transfer at hoth walls and the case of 
heat transfer at one wall.) However, the solution of the temperature 
problem (i.e., A) requires a prior knowledge of the velocity. 


TELOCITY PROBLEM 

In a manner similar to that followed by Schiller, the velocity ex- 
pression (8) is introduced into the Integrated momentum and masB equa- 
tions (5) and (7), respectively, to give the following relations between 
8, and x: 



It is convenient to introduce the dimensionless variables 

x 


16 


* _ x v _ 
a aU 

U-, 

6 * . 2 . 


e 


Re,; 




(14) 


where the equivalent diameter D 0 is equal to 4a for the duct with 
spacing 2a. Equations (12) and (13) become 
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ax*- = 

10 


3 (^1 - 1) 


TJ 


*2 


(sru* - 7) au* 


(12a) 


8 * 



(13a) 


At tlie inlet U? = 1, since TJ-, = U; ana for fully aevelopel flow 

— ■** 

= 1.5, since TT-^ = 1.5 TJ. The integration of equation (12a) is car- 
riel through numerically to = 1.5 (8* = l), at which point the velo- 
city given hy equation (8) coinciaes with the fully aevelopel parabolic 
profile. It shouia "be note! that the velocity solution given here hoes 

dIJ l 

not join smoothly to the fully levelopea parabolic profile, since — £• 

lx 

loes not approach zero as TJ* approaches 1.5. In the neighborhooa of 
the entrance section, that is, near x* = 0, the solution for the velo- 
city is 


TJ* 




(15) 


TEMPERATURE PROBLEM FOR HEAT TRANSFER AT BOTH WAHLS 


The cases of A > 8 anH A < 8 will be treatel separately. An- 
other section will be hevoted to solutions for the temperature for cases 
where the velocity profile is unchanging throughout the auct length. 

The heat-transfer parameters will be aefinea in the final section. 

Case of A > 6 . - For this case, the velocity profile in the ther- 
mal bounaary layer is flat for values of 8 ^ y < A. In the region 
0 •$ y < A, which represents the range of integration in the energy equa- 
tion (6), the velocity profile Is given by two expressions. 



0 < y ^ 5 


u 



0 ^ y ^ A 


As a consequence , the energy equation can be written as 


( 8 ) 

(16) 



(t]_ - t)u hy + TJ 




(t, - t)ay = 


r 


cJt 


y=0 


(6a) 
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vj 

y 

o 


Introduction of "the temperature expression (10) and the velocity expres- 
sion (8) yields the following dimensionless first-order differential 
equation : 



Is* I 

3 + 8 ^ ' 120 ^*3 


5 dx* 
2Pr A * 


(17) 


u V 

The quantities Uj, 8 , and x have already been defined in equation 
(14) , and 



a 


-(14a) 


The integration of equation (17) to determine the relation between A 
and x* is discussed in appendix B. Calculations were made over the 
range from A* = 0 at x* = 0 to A* = 1 (A = a) . In the neighborhood 
of the entrance, that is, near x* = 0, it was found that A* = C-^/x^, 
where C is a function of Prandtl number. 


Case of A < 8 . - Since for thiB case the velocity boundary layer 
extends beyond the thermal boundary layer, only a single expression, 
equation (8), is .required to give the velocity profile for the region 
0 4 y < A. Using the temperature expression (10) and the velocity expres- 
sion (8) in the energy equation (6) gives the following dimensionless 
differential equation: 



The integration of this equation also is discussed in appendix- B. The 
remarks already made for equation (17) concerning the range of integra- 
tion and the behavior near x* = 0 also apply to equation (18) . 


Cases of unchanging velocity profile . - The two unchanging velocity 
profiles to be considered are the parabolic and the uniform. 

Parabolic velocity profile: The equation for the parabolic profile 

with y measured from the duct wall is 



(19) 


Using this relation and the temperature expression (10) in evaluating 
the energy equation (6) leads to the following dimensionless differential 
equation : 
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J3 .*2 1 ,*3\ 3 dx* 

a (To A - 16 A J - 3? “ 

The solution for the initial condition A* = 0 at x* = 0 is 

A* 3 3 .*4 _ 3 * 

~5"-64 A “2P? X 


( 20 ) 


( 21 ) 


Uniform velocity profile: The equation for the uniform velocity 

profile is simply 

u = U 

and the solution of the energy equation (6) for initial conditions 
A* = 0 at x* = 0 is 

4 * = l/§F (22) 

Heat -transfer parameters . - The heat-transfer results will he pre- 
sented in dimensionless form utilizing the average Husselt number, which 
is defined in the usual way as 

hD e 

Hu = (23) 


For a duct with spacing 2a, the equi valent diameter D 0 is equal to 4a. 
The average heat-transfer coefficient for a length of duct x (and unit 
width) in which there is heat transfer at both walls is given by 


h = Q _ Q 
A (At) 2x(l) (At) 


(24) 


The temperature difference to be used remains to be specified. Heat- 
transfer coefficients for three definitions of At will be calculated, 
and accordingly there will be a set of Husselt numbers corresponding to 
each of the three definitions of At. These definitions are 


(At) 1 


~ (~ b w~' fc b ,x)^ 
- *1 

in — 

"W " ^b,x 


(At) 2 = (t w - t x ) 


> 


(25) 
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In the first and third definitions, t^ represents the hulk tempera- 
ture corresponding to any position x. The hulk temperature is most 
conveniently obtained from the steady-flow energy equation • 


(th z - t x ) = JL = 

0 -*■ mc„ 


Q, 


mc p p2a(l)Uc. 


(26) 


in which the kinetic energy and gravity differences have been taken to 
he negligible. 


The total heat transferred at both duct walls in a length x and 
width equal to unity is 



(27) 


The derivative in equation (27) is evaluated from the temperature ex- 
pression (10) to give 


Q 


= + 3 k( v- t 1 ) T f 


(28) 


All the quantities necessary to evaluate the IfusBelt number are 
given in equations (24) to (26) and (28); the results are as follows: 


Nu-^ = 0.25 Gz Ini 


fe) 


Uu 2 = 0.25 Gz B 


= 0.25 Gz 


B 

(1-0. 5B) 


(29a) 

(29b) 

(29c) 


The symbol Gz represents the Graetz number, which is defined by 


Be.Pr 16a^Pc_IJ 
Gz = a 




kx 


£_ _ 16 Er 


(30) 


'e *“ x 

The other dimensionless group appearing in equation (29) is given by 


B = 


2Pr 


Jo 


T* 

dx* 


(31) 


In general, B depends separately upon Prandtl number and x*. By vir- 
ture of the relation between Graetz number and x* given in equaticn 
(30), B may also be considered a function of Prandtl number and of Graetz 
number . 
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Inspection of equations (29) shows that in the most general case, 
the Husselt number will depend separately upon two other dimensionless 
groups. For example, 

"Ku = f(Gz,Pr) 

or, alternatively, 

Hu = f (x*,Pr) 

For the case of simultaneous development of the velocity and temperature 
profiles, the Husselt number depends separately on two other groups, as 
shown in equation (32) . However, for the case of the parabolic velocity 
thr oughout (and for the case of the uniform velocity profile as well), 
reference to equations (20) and (21) shows that B 1 b a function of 
Graetz number alone. Hence, for these unchanging profiles, the Husselt 
number depends on only the Graetz number. 



Although the boundary layer nature of the present analysis does hot 
allow calculations to be made for very long ducts, it is interesting to 
notice the behavior of the various Husselt numbers at great distances 
from the entrance. To this end, egression (24) for h is introduced 
into the defining equation (23) for the Husselt number to give 


Hu 


QD e 
" kA(At) 


(33) 


For very long ducts, the duct surface area A becomes very large. Also, 
the bulk temperature of the stream will approach the wall temperature, 
inspection of the definitions (25) shows that as A approaches infinity, 


At -I 


At. 


At 5 


0 

(tw-tf) 

(W 






(34) 


Equation (33), then, indicates that as A approaches infinity. Hug 
and Hu.3 approach zero. Hu]_ appears to be an indeterminate form, the 
limiting value of which has been determined in reference 3 as 7.60. 


TEMPERATURE PROBLEM FOR HEAT TRANSFER AT ORE WALL 

The temperature problem, as before, must be formulated separately 
for A > 8 and A < 8. It should be noted that this separate formulation 
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is necessary only for 5 < a;, where different formulas are required to 
describe the velocity profiles in the boundary layer and in the core. 
When 8 = a, the velocity profile is fully developed and the parabolic 
formula gives the velocity over the entire cross section. In all cases, 
the calculations are limited to A < 2a. 

Case of A > 8 . - A sketch of the development of the velocity and 
thermal boundary layers for this case is given here. 



It is seen that the thermal boundary layer, growing into the duct fhom 
the wall at temperature t^, eventually will meet the velocity boundary 
layer developing along the opposite wall. The meeting will occur when 
A = 2a - 8. In the region preceding this meeting, that is, for values 
of 6 < A '< (2a - 8), the analysis already given for the temperature 
problem with heat transfer at both walls for A > 8 is valid. 

For (2a - 8)< A < 2a, it is necessary to use three formulas to ex- 
press the velocity profile . 
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Introducing these relations together with the temperature expression 
(10) Into the energy equation (6) leads to the following result: 


l U 1 £ 3 5 + 8 A * " 120 

-s*) 2 |(i- 


.*4 


v*3 8 


,3 .* 

+ TT eA 


+ ^ f A* 2 + gA* 3 


10 


(2-6 ) (l-e) - (2 


_3_ 

4A* 


*&\ 3 


e) + | - (2-8*) 


g _ 1 f 
3 2 ^ 


( 2 - 6*) 4 


[ (1 ' a> ^ + t ^] - (2 - 5 * )5 ^ - (2 - 6 * )6 ^3 }) 


3 dx 


where 


e = 


f = 


g = 


2Pr 

A* 

4 

4 

5* 

6* 2 

4 

_ JL 


(36) 


.*2 5* 


1 ■ 


*2 


The solution of equation (36) is discussed in appendix B. 


Case of A < 6 . - The formulation of this case is identical to the 
analysis already made for the temperature problem with heat transfer at 
both walls for A < 8. 


Heat-tranBfer parameters . - The Nusselt number Nu^, based on the 
temperature difference (At)g = (t^. - t-j_) , is chosen for presentation of 
the heat-transfer results, since it is the one most conveniently used in 
calculations. If the equivalent diameter D e is taken as 4a as before, 
the expression for Nuo for the case of heat transfer at one wall becomes 
identical to equation (29b). 

PRESENTATION OF RESULTS 


Heat-Transfer Results 

Nusselt numbers for heat transfer at both walls . - The results for 
Nu^ are plotted in figure 2 as a function of Qraetz number for Prandtl 


WJ 
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numbers in the range 0.01 to 50. The curves obtained from the boundary 
layer analysis (solid lines) terminate at the Graetz numbers correspon- 
ding to A* = 1 (A = a) . The curve of reference 3 for a parabolic velo- 
city profile throughout (dot-dash line) is given for Graetz numbers be- 
low 100. Curves have been faired (dotted lines) to connect the results 
of the present analysis with those of reference 3. Figure 3 shows the 


results of Nu^ replotted with 
as parameter. 


x/D e 

Ee d 


as abscissa and Prandtl number 


The Nug results are plotted against Graetz number on figure 4 end 

2 /°e 


are replotted against 


Re; 


on figure 5. Prandtl number appears as a 


parameter on the curves. For practical applications, Hug will be the 
most convenient Nusselt number to use since it is based an the simplest 
temperature difference. The solid lines on figure 4, which represent 
the results of the boundary layer analysis, can be fitted to within 3 
percent by ‘the following equations: 


hXL 


0.644 Gz 


1/2 

Pr <p 


1+7. 



Gz , 


0.01^ Pr ^ 2 


KD, 


e 


0.664 Gz 

JF~ 


Pr 


4 11 


1+6.271 


Pr = 10, 50 


(37a) 


(37b) 


,"l/3 


The function <p appearing in equation (37a) is equal to (Pr)^' 1 ' 7 ~ for 
values of Pr 1 and is obtained from figure 11 for Pr < 1. The 
faired curves for Nu 2 shown on figure 4 are consistent with the 
faired curves drawn on figure 2 for Hu, in accordance with the 
relation *” - 


— - Gz 

Nu 2 = Tl 


4Nui 


1 - e 


Gz 


(38) 


The results for Tfu 3 are plotted in figure 6 on separate grids as a 

z/d q 

function of Graetz number and . The heat-transfer results have 

Bed 

also been presented in the foim of Uu^ because they are often given 
this way in the literature. 
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Nusselt numbers for heat transfer at one wall . - The Nug values 
for heat transfer at one wall are given on figure 7 as a function of 
Graetz number for Prandtl numbers In the range 0.01 to 50. 

Comparison of Nusselt numbers' with those of other analyses . - Two 
comparisons are made. 

(a) In reference 3, it is suggested that entrance region Nusselt 
numbers for simultaneously developing velocity and temperature profiles 
might be obtained using the results for heat transfer from a flat plate 
to a fluid having velocity U. The flat-plate results can be put into 
the following form to facilitate a comparison with the Nusselt numbers 
of the present analysis: 




0.664 Gz 


1 

2 


(39) 


For Pr ^ 1, many investigators have found that <p = Pr” ‘ . For Pr < 1, 
<p does not appear in the literature) it is calculated in appendix C of 
this report and is given in figure 11. The flat-plate results of equa- 
tion (39) are plotted in figure 8 and are compared with the Nug results 
of the present analysis for three representative values of Pr - 0.1, 2, 
and 50. 


(b) The boundary layer analysis has been used herein to calculate 
the development of the temperature profile for an unchanging parabolic 
velocity profile throughout the duct length. The Nusselt numbers (for 
heat transfer at both walls) obtained from thiB analysis are compared 
in figure 9 with those of reference 5, in which a different method of 
analysis is UBed. The values of Nu;j_ and NU3 are almost coincident 
in the range of Graetz numbers considered and appear as a single line on 
figure 9. Comparison of the Nug values from the two investigations is 
alBO made in figure 9. 


Telocity Results 

Telocity values calculated from eguation (12a) are plotted in figure 

10 as a function of = x* with the dimensionless distance from the 

aTO 

wall y/a as a parameter. Included for comparison on the same figure 
are curves representing the results of Schlichting (ref. 7). 


DISCUSSION 

Nusselt numbers for heat transfer at both walls . - Figure 2 shows 
that when Nusselt numbers for the case of simultaneously developing 
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velocity and temperature profiles are plotted against Graetz number, 
there is a separate curve corresponding to each Prandtl number. In con- 
trast , when the Eusselt numbers for the case of a parabolic velocity 
profile throughout are plotted against Graetz number, a single curve re- 
sults which is applicable for all Prandtl numbers. For any given Prandtl 
number, the Eusselt number for a flow with a developing velocity profile 
is greater than (or, at least, equal to) that for a flow where the velo- 
city profile is parabolic throughout (at a fixed Ee<j and x/D e ) because 
the velocities near the wall are higher for the developing velocity pro- 
file than for the parabolic profile. 


The curve of reference 3, based on a parabolic velocity profile 
throughout, is shown in figure 2 for the smaller Graetz numbers, where 
the boundary layer analysis can no longer be used . The previously men- 
tioned limiting value of UU]_ of 7.6 is attained in the solution of 

=4. If fully developed heat transfer is de- 


reference 3 at Gz = 


Ee^Pr 


x/K 


fined by the criterion that the limiting value of Eu^ is achieved, then 
it is seen that for a given Eeynolds number, fluids of low Prandtl number 
attain fully developed heat transfer in shorter lengths of run (i.e., 
lower x/d q ) than do fluids of higher Prandtl number. 


It is expected that at a sufficiently large distance from the be- 
ginning of the duct, the velocity conditions in the entrance region 
will cease to play an important role in determining the average Wusselt 
number. It is realized that only at infinity will the Eusselt numbers 
based on a developing velocity profile coincide precisely with those for 
a parabolic profile throughout the duct length. For practical purposes, 
however, the differences should become negligible at same finite location 
so, curves have been faired in figure 2 connecting the Eusselt numbers 
from the boundary layer analysis with those of reference 3. Although in- 
herently inexact, this extrapolation at least provides an estimate of the 
location at which the effect of velocity development is negligible. 

Figure 2 indicates that the influence of velocity development ceases at 
lower Graetz numbers for the fluids of lower Prandtl number. 


The use of the Graetz number as an independent variable allows the 
convenient comparison on one graph of results for the case of simultane- 
ous development of the thermal and velocity profiles for different 
Prandtl numbers with a Bingle curve representing the results for a para- 
bolic velocity throughout. Further, the extrapolation of the curves for 
simultaneous development of the boundary layers to join with the curve 
for the parabolic velocity 'can be made conveniently for all Prandtl num- 
bers on a single graph with Graetz number as abscissa. 

x/d 

Figure 3 is a replot of figure 2 in which — — 2. is used as abscissa 

Ee d 

and Prandtl number appears as parameter. The curves shown represent 
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the results for the case of simultaneous development of velocity and 
temperature profiles. The results for the case of a parabolic velocity 
profile throughout , which can be plotted as a single curve as a f un ction 
of Graetz number, would appear as a series of curves with Prandtl number 

as parameter if they were plotted against x ^ e 

appear on figure 3. 


Re; 


These curves do not 


The influence of Prandtl number on Nusselt number is clearly shown 

x/Dq 


in figure 3. At a given value of 


Re d 


■, the fluids of higher Prandtl 


number have higher Nusselt numbers. This result is in qualitative agree- 
ment with the usual heat -transfer correlations, which show Nusselt number 
increasing with Prandtl number for a fixed Reynolds number. The depend- 
ence of Nusselt number upon x/D 0 for a fixed Reynolds number (and a 
fixed Prandtl number) is also clearly shown in figure 3. Ab has already 
been noted, the value of x/D q corresponding to fully developed heat 
transfer decreases with decreasing Prandtl number. 


Figures 4 and 5 show trends similar to those which have been pointed 
out for figures 2 and 3. The faired curves shown on figure 4 are consist- 
ent with those drawn on figure 2. 


Nusselt numbers for heat transfer at one wall . - Nusselt numbers 
for a channel in which one wall is held at a uniform temperature while 
the other wall is insulated are shown in figure 7. The results are given 
in the form of Nu 2 , based on the temperature difference (At) 2 = t^. - tf, 
because it is the most convenient one to ubs in calculations. There 
are at present no solutions of this case for very long ducts. It is 
suggested that engineering estimates of the heat transfer for ducts 
with Graetz numbers of 3 or less be obtained from the relation 

Nu 2 = 0.5 Gz 


Comparison with other calculations . - Nusselt numbers of the present 
analysis for simultaneous development of velocity and temperature pro- 
files are compared in figure 8 with those for a single flat plate im- 
mersed in a fluid moving with a velocity TJ parallel to the plate. It 
is suggested in reference 3 that the flat-plate results might be used 
in the entrance region of a duct. For given Reynolds numbers, the agree- 
ment between the flat-plate results and those for the duct becomes 
poorer as x/D e increases. This is to be expected, since the velocity 
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of the fluid_in the core departs more and more from the free -stream 
Telocity Uoo (= U) of the plate as x/D e Increases. The agreement 
appears to he poorer for the fluids of higher Prandtl number because at 
a given Graetz number, higher values of x/D e are assocated with the 
higher Prandtl numberp at a given Reynolds number. 


A check on the boundary layer analysis was sought by calculating 
Nusselt numbers for the case of a parabolic velocity profile through- 
out the duct length. This case had already been treated in reference 3 
using Lev^que’s method for Graetz numbers above 400 and a direct numeri- 
cal integration of the energy equation for Graetz numbers below this 
value. The agreement of the curves in figure 9 is considered very good 
in view of the approximations in the boundary layer analysis and those of 
Leveque ' s method . 


Longitudinal heat flow in fluid . - The boundary layer energy equa- 
tion' (eq. (2)) is derived from the general energy equation for a constant- 
property, nondissipative flow by neglecting the longitudinal conduction 


of heat k 


S 2 t 


compared with the transverse conduction k 


S 2 i 


The 


dy 


Sx 

condition a/x « 1 is usually assumed sufficient justification for 
deleting this term. However, it is still possible that the longitudinal 
conduction can be neglected even if A/x is not small compared with 1. 
Thus, the condition A/x « 1 is sufficient, but not necessary, for de- 
leting the longitudinal conduction term. When this condition is not sat- 
isfied, the rigorous check as to whether the longitudinal conduction 
term is negligible is to obtain a solution of a problem in which it has 
not been deleted. 


For the lowest Prandtl numbers considered in this report, the con- 
dition A/x « 1 is not satisfied in the length of duct for which 
the boundary layer analysis is made. For example, inspection of figure 

x/D a 

5 shows that the solid curve for Pr = 0.01 terminates at = 0.7x10 

Ee d 


y 


at which point A = a. Since D e = 4a, this would correspond to x/A = 0.6 
for a Reynolds number of 2000. 


ax 2 


Since no solution for the energy equation containing the term 
has been obtained, there is no rigorous basis for deleting it in 


those cases where the condition — « 1 has not been satisfied. 
3 2 t 


Neglect- 


ing the term k — k for a fluid of very low Prandtl number can be justi- 
dx 2 

fied cmily on the basis of the resulting simplification of the energy 
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equation. For this reason, the results for the fluids of very low 
Prandtl number may be less accurate than the other results of this 
report . 

Telocity yrof ile . - The egression ( 8 ) chosen to describe the veloc- 
ity profile in the boundary layer had the virtue of simplicity and gave 
reasonably good results, as may be seen by comparison in figure 10 with 
the more exact values of Schlichting. It has already been noted that al- 
though the velocity solution obtained in this report coincides with the 
fully developed parabolic profile when 8 = a, its approach to the fully 
developed parabola is not asymptotic. However, it is not expected that 
any other polynomial expression for the velocity could be UBed throughout 
the entire length of the channel and asymptotically approach a parabola. 
(Basically, it would not be expected that any boundary layer analysis for 
the velocity should lead to a solution which approached the fully devel- 
oped parabola asymptotically.) 

A solution for the velocity profile which approaches the parabola 
asymptotically at infinity can be obtained by using a boundary layer 
analysis at the beginning of the duct and another method of analysis 
at great distances from the entrance. This is the method used by 
Schlichting. Such refinements in the calculation of the velocity pro- 
file would have destroyed the simplicity of the analysis for the tempera- 
ture problem which has been used in this report. 

The final check of the influence of the approximate velocity profile 
on the heat-transfer results can be made only when a more exact solution 
to the energy equation Is available. For the simpler problem of heat 
transfer from a flat plate, the author has verified that good Nusselt 
number results can be obtained with the larman-Pohlhausen method even 
though the velocity profiles are approximat'd . 

Yalldlty of thermal boundary layer calculations as A ap- 
proaches a . - Whether or not good heat-transfer results can be obtained 
from a boundary layer calculation for A values near a depends upon 
the extent of the interaction which is actually taking place in the 
duct as a result of heat transfer from tie opposite walls. It is felt 
that the boundary layer model will give good results only if this inter- 
action Is small. The values of A are in a sense arbitrary, since 
they depend to some extent upon the degree of the polynomial which is 
taken to approximate the temperature profile. For one choioe of poly- 
nomial, the condition A = a may be achieved at a position in the 
duct where the actual interaction of the heat transfer from the oppo- 
site walls is negligible. On the other hand, for another polynomial, 

A = a may occur at a location at which there is considerable inter- 
action actually taking place in the duct. 
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For a parabolic velocity profile throughout the duct length, the 
Nusselt numbers calculated from the boundary layer analysis up to A = a 
agree quite well with the result^ of reference 3. Si fact, when the 
boundary layer calculation was extended to A = 1.3a, the Nusselt numbers 
still agreed to within 3 percent. This good agreement seems to imply 
that the boundary layer model is still valid for A approaching a, at 
least when the velocity profile throughout is parabolic. 

For simultaneously developing temperature and velocity profiles, 
no other solution presently exists with which the results of the present 
analysis can be compared. So, no completely conslusive statement as to 
whether* it is valid to carry the calculation of the temperature profile 
to A = a can be made . 

Practical departures from conditions of analysis . - "Various condi- 
tions of the present analysis which may not be completely satisfied in 
practical application will be discussed. 

(a) Variable fluid properties: In some fluids, the variability of 

properties has a marked effect on the heat-transfer results. The proper- 
ty variations tend to alter the shape of the velocity profile aB well as 
to create free-convection currents. Corrections for the Nusselt number 
to take into account property changes are discussed in references 3 

and 11. 

(b) Entrance conditions: In any real apparatus, the entrance veloc- 

ity and temperature profiles will be determined by the shape of the 
transition section between the header and the duct and by the heat trans- 
fer to the fluid in the header. There may be nonzero values for the 
velocity and thermal boundary layer thicknesses at the entrance section 
associated with nonuniform profiles. 

(c) Uniformity of wall temperature: In almost any real apparatus, 

it is not expected that a uniform wall temperature can be maintained in 
the neighborhood of the entrance section; since, in addition to the pos- 
sibility of heat conduction from the duct walls to the header, the ex- 
tremely high heat -transfer coefficients near the entrance section tend 
to make the temperature nonuniform. 

(d) Heat transfer at one wall: In practice, for sufficiently long 

ducts, there may be deviations from the condition of zero heat transfer 
from one wall to the fluid even though that wall may be insulated from 
the external surroundings. This may occur when the thermal boundary 
layer growing from the heated wall contacts the opposite wall, thereby 
causing a temperature rise of the opposite wall in the region where the 
contact takes place. There will then be a conduction of heat from the 
warmer parts of the wall to the cooler parts which are located nearer 
the entrance, an a some of this conducted heat will be transferred to the 
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fluid. Except for cases of thick walls made of high conductivity material, 
this secondary heat transfer to the fluid is expected to have very little 
influence on the results of the analysis. Of course, for those ductB 
which end before A = 2a, the condition of no heat transfer to the fluid 
at one wall is achieved by insulating the wall from its external 
surroundings . 


CQNCLODIEG KEMABKS 

The simultaneous development of temperature and velocity profiles 
in the entrance region of a flat rectangular duct has been studied 
using a boundary layer analysis. The Kerman -P ohlhau sen method was UBed 
to calculate the thickness of the velocity and thermal boundary layers. 
Incompressible, constant -property fluids having Prandtl numbers in the 
range 0.01 to 50 were considered. Both duct walls were assumed to have 
the Bame uniform temperature . 

Ifusselt numbers were calculated for the length of duot for which 
the boundary layer analysis was made. It was found that when Ifusselt 
numbers were plotted against Graetz number, a separate curve for each 
Prandtl number resulted. In contrast, when the Ifusselt numbers for the 
case of a parabolic velocity profile throughout were plotted against 
Graetz number, a single curve resulted which is applicable for all 
Prandtl numbers. Eor any given Prandtl number, the Ifusselt number for 
a flow with a developing velocity profile is greater than (or, at least, 
equal to) that for a flow where the velocity profile iB parabolic through- 
out the duct length. Beyond the position where the boundary layer treat- 
ment can no longer be used, curves have been faired to connect the Ifusselt 
number results of the present analysis with those of reference 3, in 
which a parabolic velocity profile was assumed throughout the length of 
the duct. 

Calculations were also made for a duct having one wall at uniform 
temperature while the other wall was insulated . Ifusselt numbers are re- 
ported over the range of Prandtl numbers from 0.01 to 50. 

Ifusselt numbers from the present analysis were compared with those 
for a single flat plate immersed in a fluid moving with velocity 
which is equal to TJ, the average velocity in the duct. The Ifusselt 
numbers for the duct exceeded those of the plate. This difference 
increased with distance from the inlet. So that this comparison could 
be made, the IfuBBelt numbers for a flat plate for Prandtl numbers sub- 
stantially less than 1 were calculated. 

The boundary layer analysis was used to calculate Ifusselt numbers 
for a parabolic velocity throughout the duct length. These results were 
in good agreement with those of reference 3. 
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It is expected that the Nusselt numbers obtained for a flat rec- 
tangular duct will apply for an annulus for which the diameter ratio 
is near 1. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, September 30, 1954 
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APPENDIX A 
SYMBOLS 


The following symbols are used in this report: 

surface area of channel walls through which heat is Being trans- 
ferred , sq. ft 

half -spacing "between parallel walls in a duct or annulus, ft 
dimensionless group defined by eq.. (31) 
specific heat at constant pressure, Btu/(lb)(°E) 
equi valent diameter for a duct or annulus, 4a, ft 
dimensionless group defined by eq.. (36) 
dimensionless group defined By eq.. (36) 


Graetz number. 


Re d Pr 
*/ D e 


dimensionless group defined By eq.. (36) 

average heat-transfer coefficient, _^L_, Btu/(sec) (sq. ft) (°E) 
thermal conductivity of the fluid, Btu/(sec) (ft) (°F) 
mass rate of flow through cross section, lh/sec 

S) e 

average Busselt number, — 

Q \X 

Prandtl number, = _ 
k a 

static pressure, lb/sq. ft 
heat -transfer rate, Btu/sec 


diameter Reynolds number. 


LJD, 


U4a 


static temperature, °F 

bulk temperature of fluid at position x, °E 
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V 

t l 

tcD 

ir 

u i 

u* 


u„ 

u 

V 


y 

a 

A 

A* 

5 

6 * 

v 

P 

<P 


wall temperature of channel or annulus, °F 

temperature of fluid at entrance section of channel or annulus , °!F 
temperature of fluid far from surface of a flat plate, °F 


average velocity of fluid, ft /sec 
velocity in core, ft/sec 
dimensionless velocity in core, 

U 


free-stream velocity for flat plate, ft/sec 

velocity component in x-direction, ft/sec 

velocity component in y-direction, ft/sec 

coordinate measuring distance from entrance of channel, ft 

dimensionless coordinate, f SVALA 

[ a A an/ 


transverse coordinate measuring distance from channel wall, ft 

thermal diffusivity of fluid, k/c p p, ft 2 /sec 

thermal "boundary layer thickness^ ft 

dimensionless thermal boundary layer thickness. A/a 

velocity boundary layer thickness, ft 

dimensionless velocity boundary layer thickness, 6/a 

kinematic viscosity of fluid, ft 2 /sec 

density, lb/cn ft 

dimensionless group defined by equation (C4) 
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APPENDIX B 


CALCULATION PEOCEDUEE 

It was desired to use a rapid and accurate calculation procedure for 
equations (17) and (18) which also provided a continuous check. Equation 
(18) will he considered in detail. An identical procedure was used for 
equation (17). o> 

to 

If the differentiation indicated in equation (18) and the 10 

substitution 


dX* _ _ 

auT ■ 


dU* 


_1 
* ' 


(Bl) 


tl 


are made, the following first-order differential equation is obtained: 


du; 


dA 



2A* Pr U* ' 


fl A 


*2 


(B2) 


1 _ A* 2 
24 g*2 


3 

5 


A* 2 , 1 A* 5 

2 4 * c *3 




U 6 


Since 5 and U*' are simple functions of TJ* from equations (12a) 
and (13a), the numerical integration of equation (B2) to find TJ^ as a 
function of A* can be carried out without recourse to x*. This inte- 
gration was performed using the Bunge -Kutta method. In order to provide 
a continuous check at each step of the calculation, the value of 
obtained for a A* was immediately used to compute an x* from equation 
(13a) by the Bunge-Ehtta method. Alternatively, a check value of x* 
was calculated from the following integrated form of equation (18): 


x*=^ 


(5 V ‘ 24 ' f o ^(5 V " 24 6*2)“* 


(B3) 

v. 

using Simpson's rule and the values of TJ-^ and 5 obtained from the 
Integration of equation (B2). 

The starting value of dU*/dA* was obtained by using the two 
results TJ* - 1 = and A* = C , which satisfy equations (13a) 
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and (18), respectively, as x* -► 0. It should he noted that equation 
(B2) can be used only until 6* e l (u* = 1.5), at which point the 
velocity is assumed to be fully developed and U? ' becomes equal to 
zero. After the velocity profile has become fully developed and 
and S* become constants, equation (18) can be integrated directly. 

Other integration procedures were attempted for equations (17) and 
(18) involving iteration and interpolation of against x*, but these 

were found to be less accurate and no faster than the method described, 
and they did not provide a simple continuous Check.- 

It was necessary to calculate only a few points with equation (36). 
Performing the indicated differentiation would have produced a large 
number of terms requiring extensive setup time . It was decided to iterate 
equation (36), since fairly good guesses could be made for the 
values, and interpolation of P^ against x* was fairly accurate in 
the region considered. 
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APPENDIX C 


LAMINAE PQECED CONVECTION OK A HAT PLATE POE Pr < 1 

A flat plate in an infinite domain is assumed to "be alined parallel 
to a flowing stream Laying uniform, velocity and temperature upstream of 
the plate. The plate temperature t^ is held at a value different from 
that of the approaching stream. Distinct thermal and velocity boundary 
layers, A and 6, respectively, are assumed to grow along the plate. 

The case of A > 5 is considered here . A diagrammatic sketch of the 
physical model is given. 

Uco' 



The integrated momentum and energy equations for a constant-property, 
nondissipative laminar flow in a boundary layer are 


_d_ 

dx 


u(U„ - u)dy = v || 


I 

E (0. -t)*lT=aig 


y=o 


y=o 


(5a) 


(6b) 


The velocity in the boundary layer is expressed as a polynomial in 
y with coefficients that are functions of x. 


u — Ho 


KP - m 


2V5> 


0 ^y^8 


(Cl) 


This expression satisfies the following boundary conditions: 

u = cT 


d 2 u 

dy 2 


y = 0 


= 0 


u = U 

du 


7 


= 8 




= 0 
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The solution of equation (5a) can "be obtained immediately without re- 
course to the temperature "by using the velocity expression (Cl) and 
is as follows (it may also be found in ref. 12, pp. 69-70): 



(C2) 


The temperature distribution in the boundary layer is also expres- 
sible as a polynomial in y having coefficients that are functions 
of x.' 

t 

U 

This expression satisfies the following conditions: 



In obtaining the solution to equation (6b), it 1 b necessary to notice 
that in the range 0 < y ^ A, the velocity must be expressed by two 
separate functions. In the range 0 y ^ 5, expression (Cl) applies; 
while 6 < y ^ A, the velocity is given by u = Do,. Equation (6b) may 
be rewritten as 


JV 

*w 


3 y 
2 A 


y< A 


(C3) 



t)udy + Ua 


- r 

dx 


(t 0 


t)dy = a 





(C6) 


Inserting expressions (Cl) and (C3) in equation (6c) yields a first- 
order differential equation for A with a solution conveniently given 
as 


A = <p 


5 = 4.64 £ 





(C4) 
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where <p . Is the following function of Prandtl number: 


<P 


2 _ 


* ~ 8Pr " 0,4 


(C5) 


Values of (p are plotted against Prandtl number in figure 11. Also 
plotted in figure 11 is the line f 0.5 + V which is a good approx- 

\ p r 1/ y 

imation of qp for Pr ^ 0.1. 


The average Nusselt number is calculated from the definition 



Performing the indicated integration and inserting the value of h in 
equation (23a) give the final result for the average Nusselt number. 



0.664 


<p 


(C6) 
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Figure 6. - Variation of average Nusselt number Hu 3 with — and with Graetz number for flat rectangular duct. 

• Re d 

Prandtl number ap pears as parameter on curves. Fluid properties are constant. Both duct walls have same uniform 
temperature. NUj Is based on arithmetic mean temperature difference defined by equation (25). 
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Figure 7. - Plot of average Nusselt number Nug against Graetz number for flat rectangular duct with heat transfer at only 
one wall. One duot wall Is Insulated, the other Is maintained at temperature t„. Fluid properties are constant. JJu 2 
Is based on temperature difference (ty-t^) where t^ Is entering fluid temperature. 
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Figure 11. - Numerical values of function <P defined "by relation A = cp6, where A and 5 
are. respectively , thermal and velocity boundary layer thickness on single flat plate. 





